arXiv:0806.0002 

Renormalization group theory of nematic ordering 
in d-wave superconductors 

Yejin Huh and Subir Sachdev 
Department of Physics, Harvard University, Cambridge MA 02138, USA 

(Dated: May 31, 2008) 



Abstract 



We examine the quantum theory of the spontaneous breaking of lattice rotation symmetry in d- 
wave superconductors on the square lattice. This is described by a field theory of an Ising nematic 
order parameter coupled to the gapless fermionic quasiparticles. We determine the structure of the 
renormalization group to all orders in a 1/Nf expansion, where Nf is the number of fermion spin 
components. Asymptotically exact results are obtained for the quantum critical theory in which, 
as in the large Nf theory, the nematic order has a large anomalous dimension, and the fermion 
spectral functions are highly anisotropic. 
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I. INTRODUCTION 



There has been a great deal of research on the onset of a variety of competing orders 
in the hole-doped cuprate superconductors. In Ref. Q, a classification of spin-singlet order 
parameters at zero momentum was presented: such orders are able to couple efficiently to the 
gapless nodal quasiparticle excitations of a d-wave superconductor. Our focus in the present 
paper will be on one such order parameter: 'nematic' ordering in which the square lattice 
(tetragonal) symmetry of the <i-wave superconductor is spontaneously reduced to rectangular 
(orthorhombic) symmetry [2]. This transition is characterized by an Ising order parameter, 
but the quantum phase transition is not in the usual Ising universality class |l( because 
the coupling to the gapless fermionic quasiparticles changes the nature of the quantum 
critical fluctuations. Our work is motivated by recent neutron scattering observations 
of a strongly temperature (T) dependent susceptibility to nematic ordering in detwinned 
crystals of YBa 2 Cu 3 6 .45. 

An initial renormalization group (RG) analysis |l[ of nematic ordering at T = found 
runaway flow to strong-coupling in a computation based on an expansion in (3 — d) where d is 
the spatial dimensionality. More recently, Ref. Q argued that a second-order quantum phase 
transition existed in the limit of large Nf, where Nf is the number of spin components 
(the physical case corresponding to Nf = 2). We refer the reader to Ref. Q for further 
discussion on the physical importance and experimental relevance of nematic ordering in 
G?-wave superconductors. 

Here we will present a RG analysis using the framework of the 1/Nf expansion. We will 
find that a fixed point does indeed exist at order 1/Nf, describing a second-order quantum 
transition associated with the onset of long-range nematic order. However, the scaling 
properties near the fixed point have a number of unusual properties, controlled by the 
marginal flow of a "dangerously" irrelevant parameter. This parameter is v^jvp, where 
v a and vp are the velocities of the nodal fermions parallel and perpendicular to the Fermi 
surface (see Fig. [p. We will show that the fixed point has (v&/vp)* = and so the transition 
is described by an "infinite anisotropy" away from the "relativistic" fixed points found for 
other competing orders [1]. It is important to note, however, that even though the fixed point 
has (va/vf)* = 0, it is not described by an effectively one- dimensional theory of a straight 
Fermi surface; a fully two-dimensional theory is needed as is discussed in more detail in 
Section HVl The approach to this fixed point is logarithmically slow, and physical properties 
have to be computed at finite v&fvp. Our main results for the RG flow of the velocities 
are in Eqs. (13.81) and ( 13.9ft . where 6*1,2,3 are functions only of Ua/Vf which are specified 
in Eq. H3. 18[) : for v&/vp <C 1, the equations reduce to the explicit forms in Eqs. (13.221 - 
I3.24p . These (and related) equations can be integrated in a standard manner to yield the 
dependence of observables on temperature and deviation from the nematic critical point, 
and the results are in Figs. El and El 

The existence of a fixed point at (v^/vp)* = suggests that we analyze the theory 
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FIG. 1: Nodal points of the d-wave superconductor in the square lattice Brillouin zone. The 
2-component fermions are in the vicinity of the labelled nodal points and their partners at 
diagonally opposite points. 

directly in the limit of small v^/vp, without an appeal to an expansion in powers of 1/Nf. 
The nature of the small v^/vp limit is quite subtle, and has to be taken with care: it will 
be described in Section |IVl We argue there that the fluctuations of the nematic order are 
controlled by a small parameter which is ~ v&/ (NfVp), and not 1/Nf alone. Consequently, 
we believe our computations are controlled in the limit of small v^/vp even for Nf = 2. 
Indeed, the results in Eqs. (I3.22H3.2I1) . and many other related results, are expected to be 
exact as we approach the quantum critical point. 

The outline of the remainder of this paper is as follows. The field theory for the nematic 
ordering transition will be reviewed in Section [III along with a discussion of the 1/Nf expan- 
sion. The RG analysis to order 1/Nf will be presented in Section UTTI Finally, higher order 
corrections in 1/Nf, and the nature of the small v^/vp limit will be discussed in Section HVl 

II. FIELD THEORY 

We begin by review the field theory for the nematic ordering transition introduced in 
Ref . B 

The action for the field theory, S, has three components 



The first term in the action, S$, is simply that for the low energy fermionic excitations 
in the d x 2_ y 2 superconductor. We begin with the electron annihilation operator c qa at 
momentum q and spin a =T> !• We will shortly generalize the theory to one in which 
a = 1 . . . Nf, with Nf an arbitrary integer. We denote the c qa operators in the vicinity of 
the four nodal points q = (±K, ±K) by (anti-clockwise) / la , f 2a , /3a, fia- Now introduce 



S — Sq> + -\- S^)^. 



(2.1) 
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the 2-component Nambu spinors = (/i , e a bfl b ) and ^ 2a = (/2a, £abfl b ) where e ab = -e ba 
and e^i = 1 [we will follow the convention of writing out spin indices (a, b) explicitly, while 
indices in Nambu space will be implicit]. Expanding to linear order in gradients from the 
nodal points, the Bogoliubov action for the fermionic excitations of a (i-wave superconductors 
can be written as (see Fig. [T]) 



r d 2 k Nf 

= / — ^ T Y1 Yl *la (~ iUJn + v FkxT* + V A kyT X ) * la 
^ ' aj n a=l 

+ rtrf T E E $ 2« + + ^A^T 1 ) * 2a . (2.2) 

^ ' Ul n a=l 

Here u n is a Matsubara frequency, r a are Pauli matrices which act in the fermionic particle- 
hole space, k x>y measure the wavevector from the nodal points and have been rotated by 45 
degrees from q X;V co-ordinates, and Vp, v& are velocities. The sum over a in Eq. (12. 2p can 
be considered to extend over an arbitrary number Nf. 

The second term, S® describes the effective action for the Ising nematic order parameter, 
which we represent as a real field (ft. Considering only the contribution to its action generated 
by high energy electronic excitations, we obtain only the analytic terms present in the 
quantum Ising model: 



S°= d 2 xdr 



I ( a T 0) 2 + ^(V0) 2 + ^ 2 + |0 4 ]; (2.3) 

here r is imaginary time, c is a velocity, r tunes the system across the quantum critical 
point, and Uq is a quartic self-interaction. 

The final term in the action, S^^ couples the Ising nematic order, 0, to the nodal fermions, 
^ia, ^2a- This can be deduced by a symmetry analysis [lj, and the important term is a 
trilinear "Yukawa" coupling 

, N f 

= / d 2 xdr [a o J2 (*L^*la + *L^*2a) ] , (2.4) 
J 0=1 

where Ao is a coupling constant. 

We can now see that the action S describes the couplings between the <fi bosons and the 
\1/ fermions, both of which have a 'relativistic' dispersion spectrum. However, the velocity of 
'light' in their dispersion spectrums, c, Vp, v a, are not all equal. If we choose equal velocities 
with vp = f a = c, then the decoupled theory + 59 does have a relativistically invariant 
form. However, even in this case, the fermion-boson coupling Sq,^ is not relativistically 
invariant: the coupling matrix r x in Eq. (12.41) breaks Lorentz symmetry, and this feature 
will be crucial for our analysis. If we replace t x by r y in Eq. (12.41) . then we obtain an 
interacting relativistically invariant theory for the equal velocity case, which was studied in 
Ref.EI as a description of the transition from a d x 2_ y 2 to a d x 2_ y 2 + id xy superconductor. 
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A renormalization group analysis of the above action has been presented previously |l| in 
an expansion in (3 — d). No fixed point describing the onset of nematic order was found, with 
the couplings Uq and Ao flowing away to infinity. This runaway flow was closely connected 
to the non-Lorentz-invariant structure of S 1 ^ and the resulting flow of the velocities away 
from each other. In contrast, for the transition to a d x 2_ y 2 + id xy superconductor, a stable 
fixed point was found at which the velocities flowed to equal values at long scales jjj . Similar 
relativistically-invariant fixed points are also found in other cases involving the onset of spin 
or charge density wave orders at wavevectors which nest the separation between nodal points 



A. Large Nf expansion 

We will analyze S in the context of the 1/Nf expansion. Formally, this involves integrating 
out the fermions, and obtaining the resulting non-local effective action for 0. Near the 
potential quantum critical point, the non-local terms so generated are more important in 
the infrared than the terms 5? in Eq. (\2.3\\ . apart from the "mass" term r. So we can drop 
the remaining terms in S^; also for convenience we rescale — > 0/Ao and r — > NjvXq, and 
obtain the local field theory 



S = S$ + J d 2 xdr 



N 

N f r 



o=l 



(2.5) 



This local field theory will form the basis of all our RG analysis. Note that this field theory 
has only 3 parameters: r, Vp, and v&. So any RG equations can be expressed only in terms 
of these parameters, as will be presented in the following section. 

The large Nf expansion proceeds by integrating out the \& fermions, yielding an effective 
action for the nematic order 0. It is important to note that is a complicated non-local 
functional of 0, which however depends upon only r, Vp, and v&. Also, in writing out explicit 
forms for it is essential to keep subtle issues on orders of limits in mind. In particular, for 
the large Nf phase diagram, we need the effective potential for at zero k and uo. Thus we 
need an expansion of the effective potential in the regime |0|^>|A;|,|u;|-- the structure of 
this was discussed in Ref. Q, and yielded a second-order transition in the limit of large Nf. 
In contrast, for our RG analysis here, we will work in the quantum critical region, where 
(0) = 0, and so we need the functional for |0| \k\, \u\. In this case, we can expand in 
powers of to yield the formal result 



St 1 
N f 2j K 



+ 4 



r + r 2 (K))|0(A-)| 2 

I II J 5 fe K * ) V ^ K ^ K ^ ^4)0(A- 1 )0(K 2 )0(K 3 )0(K 4 ) + . . (2.6) 
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FIG. 2: Feynman graph expansion for the effective action SU. The full lines are fermion propagators, 
while the wavy lines are <f> insertions. 



Here the Ki = (ki, Ui) are 3-momenta. The functions are all given by one fermion 
loop Feynman graphs with % insertions of the external <fi vertices, as shown in Fig. [2] we 
denote the external momenta of these vertices, K iy clockwise around the fermion loop. The 
Feynman loop integrals are quite tedious to evaluate, especially for large i, and so below 
we only present explicit expressions for the needed low order terms. However, all the Tj are 
universal functions of only the momenta and v F and v A 

Our analysis will require the explicit form of T 2 (K). This can be written as 



T 2 (K) = U 2 (k x , k y , u) + U 2 (k y , k a 



(2.7) 



with the 2 terms representing the contributions of the and ty 2 fermions respectively. The 
one fermion loop diagram yields 



Yl 2 (k x ,k y ,uj) 



d 2 P rdn 



47T 



2vr 



-Tr 



r x + u) + v F (p x + k x )r z + v A (p y + k y )r x ) 1 

x t x (-iQ + v F p x r z + VAp y T x y l 
1 (uj 2 + vlkl) 



2.8) 



16v F v A (uj 2 + vpk 2 x + v\k 2 y yi 2 

Here, we have subtracted out a constant which shifts the position of the critical point. 

For T4 we will only need the cases where 2 of the four external 3-momenta vanish. The 
first is 

V A (K, -K, 0, 0) = U ia (k x , k y , uj) + U ia (k y , k x , uj) (2.9) 



where 

^-4a(k x , ky, Ui) 



Jj^J ^ Tr {r x + lu) + v F (p x + k x )r z + v A (p y + k y )r x ) x y 



X T X (-ift + V F p x T Z + V A PyT X ) 1 

1 d 2 U 2 (k x ,k y ,uj) 



2v\ 



dk 2 y 



1 (co 2 + v 2 F kl)(u 2 + v 2 F k 2 x -2v 2 A kl) 
32v F v A (u 2 + v 2 F kl + vlk 2 f/ 2 



(2.10) 
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FIG. 3: Order 1/Nf contributions to the (a) self energy Si of fermions and (b) the vertex H. 



The other case with 2 vanishing 3-momenta is 

t a (k,o,-k,o) =n 4b (k x ,k y ,u) + n 4b (k y ,k x ,u) (2.11) 

where 

/d 2 r> r dQ r f 1 2 

4^ / 2^ Tr IV ^ Q + ^ + VF ^ Px + hx)TZ + V ^ Pv + k^*) J 



(2.12) 



16v F f a (^ 2 + + v\k 2 ) 5 / 2 

From the effective action in (12. 6p . we can obtain the 1/Nf corrections to various 
observable quantities. 

For the nematic susceptibility, x<t>i we have 

i = If 2T 4 (K, -K,0,0) + T 4 (K,0,-K,0) , . 

X t> r N f J K T 2 (K)+r 1 ' ] 

Using the values in Eq. (12. 10[) and (I2.12p . we observe that the 1/Nf correction is identically 
zero. This can be traced to the arguments in Ref. Q (based upon gauge invariance and 
the fact that the coupling in Eq. ( 12. 4p is to a globally conserved fermion current) that the 
effective potential for the field <fi is unrenormalized by the low energy fermion action. Indeed, 
this argument implies that all higher order terms in 1/Nf also vanish, and that x^ 1 = r 
exactly in the present continuum field theory. So we may conclude that the susceptibility 
exponent, 7, for the nematic order parameter has the exact value 

7 = 1. (2.14) 

Note that there are 1/Nf corrections to the momentum and frequency dependence of the 
nematic susceptibility, and these will appear in our RG analysis below. 

Turning to the 1/Nf expansion for the Green's function, we write this as 

Gz (k, oj) = —iuj + vpk x r z + v^k y r x — lu), (2-15) 
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the self-energy is given by the Feynman graph in Fig. [3j 



1 f d 2 p f dil [i(p, + to) — v F (p x + k x )r z + v A (p y + k y )r x 



N f 



4tt 2 J 2tt (Q + to) 2 + v 2 F {p x + k x ) 2 + v\{p y + k y f 

1 

x 



r 2 (p,n) + r- 



(2.16) 



Finally, we will also need the 1/Nf correction to the boson-fermion vertex (the Yukawa 
couping) in Eq. (|2.4p . At zero external momenta and frequencies, this vertex is renormalized 
by Fig. [3 to 



T -| , 

Nf / 4tt 2 



d 2 P rda 



2tt 



r x (-in + v F p x r z + v A p y r x ) 1 

1 



X T 1 (-ifi + V F p x T Z + V A PyT X ) 1 T 



(2.17) 



r 2 (p, fi)+r. 

The RG equations will be obtained in the next section from Eqs. (12.131) . (12.1 6ft and (12.171) . 



III. RENORMALIZATION GROUP ANALYSIS 



We begin our discussion of the RG by describing the general structure which applies 
to all orders in the 1/Nf expansion. The RG will be based upon a consideration of the 
renormalization of the local field theory S in Eq. (12.51) . We are interested in the behavior of 
this field theory under the rescaling transformation 

k = k'e- £ 

to = Je~ l . (3.1) 

Note that we have not introduced a dynamic critical exponent z for the rescaling of the 
frequency. This is because we will allow both velocities v F and v& to flow, and this flow will 
effectively account for any anomalous dynamic scaling. Under this spacetime rescaling, we 
also have a rescaling of the fermion and boson fields 

*i ti {k,u) = ¥' li2 (fcV)e3q>Qj£ d£(4- Vf )\ 

<f>{k, to) = <f)'(k', to') exp Q J d£(5 - r] b )) . (3.2) 

Here we have allowed the anomalous dimensions rjbj to be scale-dependent, as that will 
be the case below. To implement these field rescalings, we have to determine how the field 
scales are defined. For the fermions, as is conventional, we set the scale of ^2 so that the co- 
efficients of the time derivative terms in remain unity. For the nematic order parameter, 
0, there is no "kinetic energy" term in Eq. ( 12.51) . and so we cannot use the conventional 
method. Instead, recall that the "Yukawa" coupling A was absorbed into the overall scale 
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of (j) in Eq. (12. 5p . Therefore it is natural to set the scale of <fi so that the boson-fermion 
vertex in Eq. (12.51) remains fixed at unity. This is also consonant with the fact that the 
boson "kinetic energy" comes entirely from the fermion loops. 

The RG equations now follow from the low frequency forms of the fermion self energy 
Si and the vertex H. These will have a dependence upon an ultraviolet cutoff, A. We will 
discuss below an explicit method of applying this cutoff. For now, we note that for the RG 
we only need the logarithmic derivative of the self energy and the vertex. Thus let us write 

A— £i(fc,w) = Ci(-iuj) + C 2 v F k x T z + C 3 v A k y r x 

AA S = c*-. (3.3) 

On the right hand side, we assume (as in the usual field-theoretic RG) that the limit A — > oo 
can be safely taken. Then a simple dimensional analysis shows that the Ci_4 are all universal 
dimensionless functions of the velocity ratio v^/vp. In principle, the above method can be 
generalized to any needed order in the 1/Nf expansion by taking the appropriate logarithmic 
cutoff derivative of the self energy Ei and the vertex S, as reviewed in Ref. 0- 

With the results in Eq. ( I3.3p . we can set the scaling dimensions of the fields. Considering 
the renormalization of the — iuo term in Sq,, we obtain 

Vf = -Ci. (3.4) 

Imposing the unit value of the Yukawa coupling in Eq. ( 12. 5ft we obtain 

r] b = 1 + 2C A - 2r] f = 1 + 2C 4 + 2d (3.5) 

Note the large value of r/b in the large Nf theory. This is a consequence of the fact that the 
kinetic energy is tied to the non-analytic fermion loop contribution. We can now also use 
the non- renormalization of the <ft 2 term in Eq. (12. 5p . as noted in Eq. (12.141) . to determine the 
flow equation for the "mass" r. This defines the correlation length exponent v by the RG 
equation 

% = K < 3 - 6 > 

dl v 

for small r, and we have the non-mean-field value 

u= . (3.7) 

2 - rib K J 

The results in Eq. ( 13 .3p also yield the flow equations for the velocities. By considering 
the renormalization of the fc-dependent terms in we determine that the velocity vp flows 
according 

^ = (-rj f - C 2 )v F = {C x - C 2 )v F , (3.8) 
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while the velocity v A obeys 

dv A 



M 

Clearly, the ratios of the velocity obeys 

d(v A /v F ) 



(-Vf ~ C 3 )v A = (d - C 3 )v A . (3.9) 



(C 2 -C 3 )(v A /v F ). (3.10) 



dt 

All that remains now is to determine the Ci_4 as a function of v A /v F . These functions 
also depend upon the "mass" r, but we will henceforth restrict ourselves to the vicinity of 
the critical point by setting r = 0. We reiterate that, in principle, our method allows us to 
obtain these functions to any needed order in the 1/Nf expansion. However, we will only 
obtain explicit results below to first order 1/Nf. 

The computation of C*i_4 requires the imposition of a ultraviolet momentum cutoff. We 
implement this [II] by multiplying both the Fermi and Bose propagators by a smooth cutoff 
function /C(/c 2 /A 2 ). Here JC(y) is an arbitrary function with /C(0) = 1 and which falls off 
rapidly with y at y ~ 1 e.g. JC(y) = e~ y . We will show below that the RG equations are 
independent of the particular choice of JC(y). 

Let us now discuss the evaluation of C\- 3 . Schematically, we can write for the fermion 
self energy 

Sl(A ' ) = /0 F(p + A ' )G(P KTMfi (3n) 

where K = (k, uj) and P = (p, fl) are 3-momenta, and F and G are functions which can be 
obtained by comparing this expression to Eq. (12.161) . Also notice that (at r = 0) F and G 
are both homogenous functions of 3-momenta of degree -1. Expanding to first order in 
(fjL is a spacetime index) we have 

- *>J B? G ^ © + F ^ © ©] « 312 » 

where p M is the three vector (0,p x ,p y ) (while P M = (Q,p x ,p y )). Now taking the A derivative 



(3.13) 



Now we convert to cylindrical co-ordinates in spacetime by writing P M = yA(vpx, cos 9, sin i 
and integrate over y, x, and 9. Also, let us define P M = (vpx, cos 9, sin 9) and p^ 
(0, cos 9, sin 6). We use the homogeneity properties of the functions F and G to obtain 



POO 

- Ap,F{P)G{P) / y 3 dy {K {y 2 ) K" (y 2 ) + K! 2 (y 2 ) } 
Jo 



(3.14) 
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Now the integrals over y can be evaluated exactly by integration by parts, and the final 
result is independent of the precise form of K{y): 

In a similar manner, we can write the result for the vertex 5 in Eq. (12.171) in the schematic 
form 

s ^+/0w(£>' (3i6) 

where H(P) is a homogenous function of P of degree -3. Proceeding as above we obtain the 
analog of Eq. (IXT5]) 

d „ v F ro ° r27T 



A iK E= ^L dx J dBH(p) - (317) 

We can now use the above expressions to obtain explicit results for the constants Ci_ 4 
at order 1/N f . For Ci_ 3 we combine Eqs. ( l2~T6l) . (1331) . (l3TTj) . and fl3~T5|) to obtain 

r _ 2(v±/v F ) [ 2 * ft (x 2 -cos 2 g-(W^) 2 sin 2 g) 

1_ tt»JV> 7_ 0O ^ 7o ^(x 2 + cos 2 ^ + (W^) 2 sin 2 ^^ lX ' yj 

Ti Nf > /_ 00 ,/ (x 2 + cos 2 + (wa/^f) 2 sin 2 #) 2 

7T 6 Nf ,/_ 00 y (x 2 + cos i 6 l + (v A /v F ysm z 9y 

where 

- i , „ s x 2 + sin 2 x 2 + cos 2 6 

g-\x, 9) = ===== + ===== (3.19) 

a/x 2 + (v&/v F ) 2 cos 2 6* + sin a/x 2 + cos 2 9 + (v&/v F ) 2 sin 6 1 

is the propagator inverse. All the integrals in Eqs. ( 13. 18ft are well-defined and convergent 
for all non-zero va/v f , and have to be evaluated numerically. It is also evident that the 
results are functions only of v&/v F . The results of a numerical evaluation are shown in 
Fig.H 

Similarly, we can combine Eqs. (12.171) . (13.31) . (I3.16p . and (13.171) to deduce that 

C 4 = -C 3 , (3.20) 

at leading order in the 1/Nf expansion. 

With the values in Eq. ( 13.181) . we are now in a position to numerically solve the RG flow 
for the velocity ratio v&/v F in Eq. (|3.10|) . We plot the right-hand-side of this flow equation 
in Fig. [51 Note that the only zero of the flow equation is at Va/v f = 0, and that this 
fixed point is attractive in the infrared. So for all starting values of the parameter v&/v F , 
the flow is towards v&/v F — > for large t. At large enough scales, it is always possible 
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FIG. 4: Plots of the functions Ci_ 3 in Eq. (|37l8l) . 
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FIG. 5: Right-hand-side of the RG flow equation Eq. (|3,10p for v A /v F . 



to approximate the constants C*i_3 in Eqs. (13.181) by their limiting values for small v A /v F . 
Evaluating the integrals in Eq. (13.181) in this limit, we obtain 

Cl = -0.4627^M + o ((v A /v F f) 
JSf 



C, -0.3!79^M + O(K/^) 3 ) 



C, 



TV 



2 \n(v F /v A ) - 0.9601) + o {(v A /v F ) 3 ) 



(3.21) 



The logarithmic dependence on v A /v F arises from a singularity in the integrand at x = 
and 9 = ±7r/2: this corresponds to an integral across the underlying Fermi surface over uo 
and k x where the inverse fermion propagator ~ —iuo + v F k x r z . Inserting these values into 
the flow equations (13.81) . (13.91) . and 1 13. lOft we obtain an explicit form of the RG flow for 
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small va/vf- 



~d£ 

dv F 

d(v A /v F ) 



7T Z 



\n(v F /v A ) - 0.4974 



NfVp 



-0.1148^ 

N f 



— In (v F /v A ) 



(3.22) 
(3.23) 
(3.24) 



0.6122") MVF) \ 

The right-hand-side of Eq. (13.241) has a zero for va/v f ~ 2; this is spurious as this equation 
is valid only for v A /v F <C 1. The full expression for the right-hand-side, valid for arbitrary 
v A /v F , is plotted in Fig. and this makes it clear that the only zero of the beta function is 
at v A /v F = 0. 

The results of a numerical integration of Eqs. (13.81) . (13. 9p and (I3.10p are shown in Figs. El 
EJ and El We start the RG integration at £ = corresponding to a temperature T = T 
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FIG. 6: Integration of Eqs. ([3JiD, (USD an d MP from <= to £ for N f 
velocities v F and v^. Results above are for v^/v F = 1. 



2, starting from the 



at the nematic critical point, at which point the velocities have the values v° F and v\. 
Integrating to the scale £ yields the velocities v F (£) and v A {£) at a temperature T = T e~ e . 
The results depend upon the initial value of the velocity ratio, v A /v F , and are shown in the 
figures for v A /v F = 1, 0.1, 0.05. These results can also be converted into dependence on the 
deviation from the nematic critical point, r, by integrating Eq. (13. 6p . 

In the asymptotic low temperature regime {£ — * oo), we have v A /v F — > 0, and so we can 
integrate the asymptotic expression in Eq. (13.241) . This yields 

^ 1 (3.25) 



VA 
V F 



8 £\n[0. 3809 £/N f ] 
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FIG. 7: As in Fig. El but for v° A /v° F = 0.1. 
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FIG. 8: As in Fig. ©but for v° A /v° F = 0.05 



The T dependence follows by using £ = \n(T /T). Using the result for va/vf in Eq. ( 13. 25ft 
and integrating Eq. (13.231) . we find that vp only has a finite renormalization from its bare 
value as t — > oo. So the decrease in v^jvp as £ — > oo (i.e. as T — > 0) comes primarily from 
the decrease in t> a- This is also clear from the plots in Figs. [ME1 
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We can also apply these results to obtain the T = evolution of the velocities as a 
function of the tuning parameter, r, across the nematic ordering transition. In this case 
we use Eq. (I3.6p . along with the fact that v — > 1 as i>a/Vf — > 0. Then we deduce that 
I = ln(l/|r|) to leading logarithm accuracy, and so obtain the r dependence of the velocities 
from Figs. [MHl When both T and r are non-zero, we can use I = ln(max(|r|, T)) to leading 
logarithmic accuracy. Note that this predicts a minimum in v^/vp as r is tuned across the 
nematic ordering transition, with the minimum value of order l/ln(l/T). 



IV. THEORY FOR SMALL v A /v F 

The RG analysis in Section HTT1 has shown that the (3 functions at order 1/Nf are such 
that the flow is towards the strong anisotropy limit, va/vf — > at large scales. Here we 
will show that this conclusion holds at all orders in the 1/Nf expansion. Indeed, for small 
va/vf we will show that va/vp can itself be as the control parameter for the computation, 
and Nf is not required to be large. In other words, the RG flow equations in Eqs. (13.221 
I3.24p are asymptotically exact, even for the physical value of Nf = 2. All corrections to 
Eq. (I3.22H3.24"|) are higher order in v/^/vp- 

Our conclusions follow from an examination of the structure of the action for the 
nematic order in Eq. (12.61) in the limit of i>a/Vf — > 0. As we will shortly verify, the natural 
scale for the fluctuations of are at frequency and momenta with u ~ vpk x ~ Vpk y . Under 
these conditions, we can just take the v a — > limit of the expressions in Eqs. (12.81) . (12.101) . 
and (12.121) . This limit exists, and to quadratic order in we have the following effective 
action as Va — > (at criticality, r = 0): 



l§VpV A J 4tt 2 J 2tt v ' ; 



u 2 + v%kl + xjij 1 + Vpk 2 . 



+ ... (4.1) 



Rescaling momenta k — > kjvp, we notice that the prefactor in Eq. (14.11) is proportional to 
the dimensionless number N^vp/v^- Thus each propagator of <fi will appear with the small 
prefactor of v&/(NfVp). Further, as claimed above, each propogator has the typical scales 
to ~ v F k x ~ vpky. 

These considerations are easily extended to terms to all orders in in Eq. ( 12.61) . The 
full expression for involves the sum of the logarithms of the determinants of the Dirac 
operators of fermions in a background <fi field. The expansion of these determinants in 
powers of involves only terms with a single fermion loop. For each \&i loop we rescale the 
fermion momentum p x — > p x /vf and p y — > Py/v^. Similarly for each \I>2 loop we rescale 
the fermion momentum p x — > p x jv^ and p y — > p y /vp. In both cases, the integral over the 
fermion loop momentum yields a prefactor of I/^f^a)- For the vertices emerging from 
the loops, we have momenta uo ~ vpk x ~ vpk y , as noted above; when this external 
momenta enters the fermion loop, we have k x ~ p x , the typical x-momentum. However, 
k y p y ~ 1/fA) the typical y-momentum. Consequently, we may safely set k y = 
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within the fermion loop, and there is no further dependence upon va in the loop integral. 
Similarly, in the ^2 loops, we may set k x = within the fermion loop. To summarize, the 
typical <p momenta are k x ~ l/vp, k y ~ 1/Vf, the typical momenta are p x ~ 1/vp, 
p y ~ 1/fAj and typical ^2 momenta are p x ~ 1/^a, Pj, ~ 1/Vf- Further, the resulting 
effective action for has the form 

Sff, = S[(j)(k, to); Vpk x , Vpky, to] + terms higher order in v^/vp- (4.2) 

VfVa 

Here S is a functional of <p which is independent of v a, whose expansion in has co-efficients 
which depend only upon VFk x , vpk y , and to. It can be verified that Eq. (14. II) and the quartic 
terms in Eqs. (12 . 10[) and (12.121) are indeed of the form in Eq. (I4.2p . 

From Eq. (14.21) . it is evident (after rescaling momenta k — > k/vp and momentum space 
fields — > v F (p) that the natural expansion parameter controlling fluctuations is v a/ (Nvf). 
This is the main result of this section. 

It remains to understand the \to.(vf/va) factor in Eq. (13.241) . If we compute observable 
<p correlations under the reduced effective action S, the resulting expansion in powers of 
v&/(Nvf) leads to Feynman integrals which are not necessarily infrared or ultraviolet finite. 
Because the action S is free of any dimensionful coupling constant, the Feynman graph 
divergences can at most be powers of logarithms. These divergences are cutoff by using 
the full fermion propagators, including the terms proportional to v a times the momenta 
of the boson propagators. This leads to factors of ln(w f / 'v a) , as was the case in obtaining 
Eq. (l3T24j) from Eq. (13481) . 



V. CONCLUSIONS 

This paper has described the RG properties of the field theory jl] in Eq. f 1 2 . 5 1) for nematic 
ordering. It was previously noted [4| that the effective potential for the nematic order, 
0, remained unrenormalized upon integrating out the nodal fermions ^1,2- This happens 
because the <fi couples to a conserved fermion current, and a spacetime- independent <fi can 
be 'gauged away' applying a gauge transformation to the fermions. One consequence of 
this non-renormalization is that the susceptibility exponent, 7, takes the simple value in 
Eq. (I2.14p . However, non-trivial renormalizations of the field scale of both <p and \&i 2 are 
still possible, along with renormalizations of the velocities, and we have shown here that 
this leads to an interesting RG flow structure, with non-mean-field exponents. 

It is interesting to note here the similarities to the RG flow of supersymmetric field theories 
[sj]. There the effective potential also remains unrenormalized (albeit for very different 
reasons), but the 'wavefunction renormalizations' lead to many non-trivial RG fixed points. 

Our main result was that the transition is described by a fixed point in which the fermion 
velocities at the nodal points have a ratio which approaches a fixed point with va/vf — > 
logarithmically slowly. However, it is not valid to set the pairing-induced velocity = in 
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the computation, and so deal with a metallic Fermi surface. For the case where the fermion 
dispersion is ~ (v F k 2 + v A kyY^ 2 , the typical fermion momenta contributing to the critical 
theory scale as k x ~ 1/vp and k y ~ 1/fA) and so the full form of the Bogoliubov quasi- 
particle dispersion is important. The flow of the velocities is described by Eqs. (I3.22H3.24|) 
for va/vf <C 1, and these equations are believed to be asymptotically exact. 

Unlike the fermions, the fluctuations of the nematic order, 0, are isotropic in momenta. 
These are described asymptotically exactly by Eq. (14. ip . Note that the propagator is very 
different from free field, and has a large anomalous dimension, as was found in the Nf = oo 
theory [4 Experimental detection of this unusual spectrum in e.g. inelastic X-ray scattering 
would be most interesting. 
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